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We study phenomenologically the scenario in which the scalar top quark is 
lighter than any other standard supersymmetric partner and also lighter than 
the top quark, so that it decays to the gravitino via t — > W + bG. In this case, 
scalar top quark events would seem to be very difficult to separate from top 
quark pair production. However, we show that, even at a hadron collider, it is 
possible to distinguish these two reactions. We show also that the longitudinal 
polarization of the final W + gives insight into the scalar top and wino/Higgsino 
mixing parameters. 
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1 Introduction 



Supersymmetry has been studied for a long time as the possible framework for elementary 
particle theories beyond the standard model JE], |J. It provides a natural solution to the 
hierarchy problem, allowing a small value, in fundamental terms, for the weak interaction 
scale. It also allows the measured values of the standard model coupling constants to be con- 
sistent with grand unification. Still, if Nature is supersymmetric, some new interaction must 
spontaneously break supersymmetry and transmit this information to the supersymmetric 
partners of the standard model particles. Two different approaches have been followed to 
model supersymmetry breaking. The first is the idea that supersymmetry breaking is trans- 
mitted by gravity and supergravity interactions f§]. In these scenarios, the supersymmetry 
breaking scale a/F is of the order of 10 11 GeV. This large value implies that gravitino in- 
teractions are extremely weak, and that the gravitino has a mass of the same size as the 
other supersymmetric partners. In this class of models, the lightest supersymmetric parti- 
cle (LSP), which is the endpoint of all superpartner decays, is most often taken to be the 
superpartner of the photon, or, more generally, a neutralino. 

The second approach uses the gauge interactions to transmit the information of supersym- 
metry breaking to the standard model partners H [| [7|- In these gauge- mediated scenarios, 
the supersymmetry-breaking scale \fF is typically much smaller than in the gravity-mediated 
case, so that the gravitino G is almost always the LSP. All other superpartners are unstable 
with respect to decay to the gravitino, though sometimes with a lifetime long on the time 
scale relevant to collider physics. 

In gauge-mediated scenarios, direct decay to the gravitino is hindered by a factor 1/F 
in the rate. Thus, attention shifts to those particles which have no allowed decays except 
through this hindered mode. Such a particle is called a next-to-lightest supersymmetric 
particle (NLSP). Any of the typically light superpartners can play the role of the NLSP, and 
the collider phenomenology of a given model depends on which is chosen. For example, if the 
gaugino-like lightest neutralino x° is the NLSP and decays inside the collider, supersymmetry 
reactions will end with the decay x° ~ > jG, producing a direct photon plus missing energy. 
Other common choices for the NLSP are the lepton partners and the Higgs boson. More 
involved scenarios are also possible [|J. 

In this paper, we consider the possibility that the the lightest scalar top quark (stop, or 
ti) is the NLSP of a gauge-mediation scenario ||. It is typical in supersymmetric models that 
the stop receives negative radiative corrections to its mass through its coupling to the Higgs 
sector. In addition, the mixing between the partners of the ti and tn is typically sizable, 
and this drives down the the lower mass eigenvalue. It is not uncommon in models that the 
lighter stop is lighter than the top quark, and it is possible to arrange that it is also lighter 



than the sleptons and charginos JED], [TTJ. The existence of this possibility, though, poses 
a troubling question for experimenters. In this scenario, the dominant decay of the lighter 
stop would be the three-body decay t — > bW + G. The G is not observable, and the rest of 
the reaction is extremely similar to the standard top decay t — > W + b. The cross section for 
stop pair production is smaller than that for top pair production at the same mass. Thus, it 
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is possible that the top quark events discovered at the Tevatron collider contain stop events 
as well. How could we ever know? In this paper, we address that question. 

Our strategy will be to systematically analyze the three-body stop decay. This decay 
process is rather complex, since the G can be radiated from the partners of t, b, or W, 
and since both the top and the W partners can be a mixture of weak eigenstates. For the 
application to the Tevatron, one must take into account that the center-of-mass energy of 
the production is unknown, and that the detectors can measure only a subset of the possible 
observables. Nevertheless, we will show that two observables available at the Tevatron can 
cleanly distinguish between top and stop events. The first of these is the mass distribution 
of the observed b jet plus lepton system which results from a leptonic W decay. The second 
is the W longitudinal polarization. We will show that the first of these observables gives a 
reasonably model-independent signature of stop production, while the second is wildly model- 
dependent and can be used to gain insight into the underlying supersymmetry parameters. 

This paper is organized as follows: In Section 2, we set up our basic formalism and state 
our assumptions. In Section 3, we analyze the stop decay rate and the bW and bi mass 
distributions. In Section 4, we present the W longitudinal polarization in various models. 
Section 5 gives our conclusions. 

2 Formalism and assumptions 

In this section, we define our notation and set out the assumptions we will use in analyzing 
the stop decay process. Our calculation will be done within the framework of the minimal 
super symmetric standard model (MSSM) with R-parity conservation. We will not consider 
any exotic particle other than those required in the MSSM. 

Our central assumption will be that the lighter stop mass eigenstate t\ is lighter than the 
top quark and also lighter than the charginos and the b superpartners, while the gravitino is 
very light, as in gauge- mediation scenarios. Under these assumptions, what would otherwise 
be the dominant decay ti — > tG is forbidden kinematically, so that the dominant stop decay 
must proceed either by t\ — > bW + G or by t\ — > cG. In the MSSM without additional flavor 
violation, quark mixing angles suppress the decay to c by a factor 10~ 6 . That suppression 
makes this decay unimportant except near the boundary of phase space where m ~ raw- 
For this reason, we will ignore that decay in the rest of the paper. 

If the mass of the t\ were larger than the mass of the top quark, the t\ would decay 
entirely through t\ — > tG. All observable characteristics of this decay would be exactly 
those of top quark pair production, except that the two emitted gravitinos would lead to 
a small additional transverse boost. For such a heavy stop, the production cross section is 
less than 10% of that for top quark pair production. Nevertheless, this process might be 
recognized from the fact that the top quark and antiquark would be given a small preferential 
polarization, for example, in the t^L helicity states if the t\ is dominantly the partner of 
t R . The methodology of the top polarization measurement has been discussed in detail in 
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the literature [JT2 , so we will not analyze this case further here. 



To analyze the case in which t\ is lighter than the top quark, we begin by considering 
the form of the scalar top quark mass matrix. Including the effects of soft breaking masses, 
Yukawa couplings, trilinear scalar couplings, and D terms, this matrix can be written in the 
tji, t L basis as 

M? = ( m L m t (A t + fj, cot (3) \ 

1 \m t {A t + fMcotp) m\ L ) ' K ' 

where A t , fi, m t , and tanf3 denote, respectively, the trilinear coupling of Higgs scalars and 
sfermions, the supersymmetric Higgs mass term, the top quark mass, and the ratio of the 
two Higgs vacuum expectation values. The masses m 2 R and arise from the soft breaking, 
the D term contribution, and the top Yukawa coupling as follows: 

2 

m~ = mfr + m\ H — sin 2 6 w m 2 z cos 2f3 

,; . ., t ., ! .; H .,.. ) . i ., ! , 

where 6 W denotes the weak mixing angle and rriz is the Z° boson mass. The soft breaking 
masses m 2 ^ and m 2 ^ are more model-dependent. In many models, these masses are derived 
from flavor-blind mass contributions by adding the effects of radiative corrections due to the 
top-Higgs Yukawa coupling Af. These corrections have the form 

m l 3 ~ m l- 2 >?J , m % 3 ~ m | - )?Ji ( 3 ) 

where the function / denotes a one-loop integral. The extra factor 2 in the expression for the 
is due to the fact that loop diagram contains the Q and Higgs isodoublets. From this 
effect, we expect that m% < m% . One should note that there is a flavor- universal positive 
mass correction due to diagrams with a gluino which combats the negative correction in (|3|). 

The lightest stop mass eigenstate t\ and its mass m 2 are easily obtained by diagonalizing 
the stop mass matrix (|TJ). One finds 

t\ = cos 9ttL + sin 9ft r 
t 2 = sin9 t t L - cos 9 t t R 

m 2 = ^{m 2 R + m\ - yj {m\ h - mfj 2 + Am 2 t {A t + Jl cot f3) 2 } 

m t (A t + /icotp) 

tarUf t = — 2 ^ — • ( 4 ) 

(mf R - m{) 

In these formulae, 9$ denotes the stop mixing angle and is chosen to be in the range —tt/2 < 
9i < 7r/2. The relations @ demonstrate the two mechanims mentioned in the introduction 
for obtaining a small value of mj x : First, the radiative correction (Q) could be large due to 
the large value of A*; second, the left-right mixing could be large due a large value of A t . 
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From here on, however, we will take m and 6t to be phenomeonological parameters to be 
determined by experiment. 

Since the final state of the three-body ti decay includes the W + , our analysis must include 
the super symmetric partners of W + and H + , the charginos. In the MSSM, these states are 
mixtures of the winos and the Higgsinos h . In two-component fermion notation, the 
left-handed chargino fields are written 

Cf = (w + , ih+), Cf = (w-,ih-). (5) 

In this basis, the chargino mass matrix is 

\— Vzmw cos p [A / 

where is the soft breaking mass of the SU(2) gaugino, and fi is the supersymmetric 
Higgs mass. The matrix M + is diagonalized by writing M + = (V_) T DV+, where V+, V- are 
unitary; then the mass eigenstates are given by 

Xt = V +ij C; , xT = V-ifiJ . (7) 

To be consistent with the assumption that the t\ is the NLSP, we will consider only sets of 
parameters for which the mass of the t\ is lower than either of the eigenvalues of M + . 

We analyze the couplings of superparticles to the gravitino by using the supersymmetry 
analogue of Goldstone boson equivalence. The gravitino obtains mass through the Higgs 
mechanism, by combining with the Goldstone fermion (Goldstino) associated with sponta- 
neous supersymmetry breaking. When the gravitino is emitted with an energy high compared 
to its mass, the helicity h = ±| states come dominantly from the gravity multiplet and are 
produced with gravitational strength, while the h = ±| states come dominantly from the 
Goldstino. In the scenario that we are studying, the mass of the gravitino is on the scale of 
keV, while the energy with which the gravitino is emitted is on the scale of GeV. Thus, it 
is a very good approximation to ignore the gravitational component and consider the grav- 
itino purely as a spin | Goldstino. From here on, we will use the symbol G to denote the 
Goldstino. 

The coupling of one Goldstino to matter is given by the coupling to the supercurrent |L3| 

6C = ~72F 9fldcJ " + 72F J ^ CdflG * ' (8) 

where \/F is the scale of supersymmetry breaking and c = — ia 2 . The supercurrent takes 
the form 



r = y/2(J^D v <f> m i/)-V2i(^J 



-ga"al)*\*(j) - ia Xa F Xa a»c\* , (9) 
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summed over all chiral supermultiplets (<f>,if)) and all gauge supermultiplets (A^, A). In this 
equation, W is the superpotential and g is the gauge coupling. All of the various terms in 
this equation actually enter the amplitude for the three-body stop decay. 

It is a formidable task to present the complete dependence of the properties of the three- 
body stop decay on the various supersymmetry parameters. We will present results in this 
paper for the following four scenarios, which illustrate the range of possibilities for the wino- 
Higgsino mixing problem: 

1. a scenario in which the lightest chargino is light and wino-like: m 2 = 200 GeV, \x = 1000 
GeV, 

2. a scenario in which the lightest chargino is light and Higgsino-like: m 2 = 1000 GeV, 
fi = 200 GeV, 

3. a scenario in which the lightest chargino is light and mixed: m 2 = \x = 260 GeV, 

4. a scenario in which the lightest chargino is heavy: m 2 = /x = 500 GeV. 

Within each scenario, we will vary other parameters such as m, sin 6^, and tan/? in order to 
gain a more complete picture of the t\ decay. 



3 Characteristics of the stop decay 

Using the Goldstino interactions from (|8|) and the gauge and Yukawa interactions of the 
MSSM, we can construct the Feynman diagrams for t\ — > bW + G shown in Figure [1|. These 
diagrams include processes with intermediate t, xf, and b particles, plus a contact interaction 
present in (||). 

It is useful to think about building up the complete amplitude for the stop decay by 
successively considering a number of limiting cases. In Figure [TJ, we have drawn the diagrams 
using a basis of weak interaction eigenstates. 

The first property to be derived from these amplitudes is the stop decay rate. It is always 
an issue when an NLSP decays to the gravitino whether the decay is prompt on the times 
scales of particle physics, or whether the NLSP travels a measureable distance from the 
production vertex before decaying. Taking into account the 3-body phase space and the fact 
that the amplitude is proportional to 1/F, we might roughly estimate the decay amplitude 
as 

~ a w {m-m w ) 7 
T{tl) ~ 1028ttX^ ' (10) 

where a w = g\j^ is the weak-interaction coupling constant. By this estimate, a value of 
\fF smaller than 100 TeV would give a prompt decay, with cr < 1 cm. 

In Figure ^| we show the result of a complete calculation of the decay rate in the four 
scenarios listed at the end of Section 2. In all four cases, we have chosen the parameter values 
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Figure 1: Feynman diagrams for the process t\ — > bWG. All diagrams are drawn in terms 
of 2-component fermion notation. G denotes the Goldstino. The label on the W/h internal 
lines labels the vertex with which the chargino couples to the top quark. 
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Figure 2: Decay length for the lightest scalar top quark t\ as a function of its mass in four 
different scenarios. The parameters \fF = 30 TeV, tan/5 = 1.0, = 300 GeV, and 
sin^t = —0.8 are assumed in all four scenarios. 
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\/F = 30 TeV, tan/3 = 1.0, m~ bL = 300 GeV, and sm6 t = —0.8. The complete calculation 
reproduces the steep dependence on the stop mass m which is present in (|T0|) . and shows 
that the normalization is roughly correct. Since V varies as the fourth power of y/F, one 
can arrange for a short decay length by making \/F sufficiently low. For m ~ 160 GeV, the 
choice \/F = 30 leads to a decay length cr of about 1 cm. We have found that the decay 
length is quite insensitive to all of the other relevant parameters. The variation between 
scenarios or within a given scenario is less than a factor of 2. From here on, we will analyze 
the t± decay as if it were prompt. But it is clear from the figure that, if \fF is as low as 30 
TeV, stop decays will be identifiable by their displaced vertices in addition to the kinematic 
signatures discussed in this paper. 

The final state of the three-body stop decay is essentially the same as ordinary top decay, 
since the stop produces a b jet, a W boson, and an unobservable G. How, then, can we 
distinguish the tt and tit\ production processes? The most straightforward way to approach 
this problem is to analyze the observable mass distributions of t and t\ decay products. If we 
could completely reconstruct the W boson, the invariant mass of the bW system would peak 
sharply at m t in the case of t decay, and would have a more extended distribution below 
the stop mass m in the case of t% decay. However, in the observation of top events at the 
Tevatron, the analysis cannot be so clean. Events from tt production are typically observed 
in the final state in which one W decays hadronically and the second decays to Iv. Then 
the final state contains an unobserved neutrino. If there is only this one missing particle, 
the event can be reconstructed. But the events with t\ contain two more missing particles, 
the Gs, which potentially confuse the analysis. 

Fortunately, it is possible to discriminate t from t% events by studying the invariant mass 
distribution of the directly observable b and lepton decay products. For top decays, the 
distribution in the 6-lepton invariant mass m(eb) (quoted, for simplicity, for m;, = 0) takes 
the form 

1 dT 12m(efe) m 2 

F^py = 2(1 - m 2 w /m 2 )(2 + m 2 w /m 2 ) ~ ~ ^ + T t^ v V) ' (U) 

where y = m 2 (eb)/(m 2 — m^). As is shown in Figure || this distribution extends from 
m(eb) = mb to a kinematic endpoint at m(eb) = 155 GeV, and peaks toward its high end, at 
about m(eb) = 120 GeV. On the other hand, in t\ decay, not only does the m(eb) distribution 
have a lower endpoint value, reflecting the value of m < m t , but it also peaks toward the 
low end of its range. Figure ^| shows two typical distributions of m(eb), corresponding to 
stop masses of 130 and 170 GeV. The corresponding distributions of the b-W invariant mass 
m(bW) are also shown for comparison. 

A remarkable feature of Figure |^ is that the m(eb) distributions from top and stop 
decay remain distinctly different even in the limit in which the stop mass m approaches m t . 
Naively, one might imagine that the stop decay diagrams with top quark poles, (a) and (b) 
in Figure [I], would dominate in this limit and cause the stop decay to resemble top decay. 
Instead, we find that the top pole diagrams have no special importance in this limit. If Eq 
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Figure 3: Typical eb and Wb invariant mass distributions for t\ decay. The distributions are 
shown at two different assumed t\ masses, = f30 (dotted lines) and f 70 GeV (dashed 
lines), under the scenario (f). The other parameters are chosen as tan/3 = 1.0, m~ hh = 
300 GeV, and sin 6*t = —0.8. The solid line shows, for comparison, the m e & spectrum for the 
standard top quark decay. 
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is the G energy, the top quark pole gives an energy denominator 1/Eq, but this is cancelled 
by a G emission vertex proportional to (Eq) 3 ^ 2 . 

In Figure |], we show the variation of the the distribution of m(eb) and m(bW) according 
to the choice of the supersymmetry parameters. The five curves in each group correspond 
specific parameter choices in the four scenarios listed at the end of Section 2, plus an addi- 
tional choice in scenario (3) corresponding to the case of a pure tjj (9 t = 0). The distributions 
for a given value of m are remarkably similar. Presumably, the shape of these distributions 
is determined more by general kinematic constraints than by the details of the decay ampli- 
tudes. The only exception to this rule that we have found comes in the case where the t\ is 
dominantly ti and the w exchange process is especially important. 

From these results, we believe that the t\ production process can be identified by measur- 
ing the distribution of m(eb) in events that pass the top quark selection criteria. The mass 
of the ti can be estimated from this distribution to about 5 GeV without further knowledge 
of the other supersymmetry parameters. 



(12) 



4 Longitudinal W polarization 

One of the characteristic predictions of the standard model for top decay is that the final- 
state W bosons should be highly longitudinally polarized. Define the degree of longitudinal 
polarization by 

_ r(w ) 
r r(aii) ' 

Then the leading-order prediction for this polarization in top decay is 

n = 1 , 9 « 0.71 . (13) 

l + 2m 2 w /m 2 1 ; 

We have seen already that the configuration of the final bW + system in stop decay is quite 
different from that in top decay. Thus, it would seem likely that the longitudinal W polar- 
ization would also deviate from the characteristic values for top. We will show that the value 
of r in stop decay typically differs significantly from (pT3[), in a manner that gives information 
about the underlying supersymmetry model. 

The measurement of the polarization r at the Tevatron has been studied using the tech- 
nique of reconstructing the W decay angle in single-lepton events from the lepton and neu- 
trino four- vectors 0, 15| . An accuracy of ±0.03 should be achieved in the upcoming Run 



II. This technique, however, cannot be used for stop events, since the missing momentum 
includes the G"s as well as the neutrino. However, one can also measure the longitudinal W 
polarization from the W decay angle determined by using the four-vectors of the two jets 
assigned to the hadronic W in the event reconstruction. It is not necessary to distinguish 
the quark from the antiquark to determine the degree of longitudinal polarization. 

What value of r should be found for light stop pair production? In Figures | and |], we 
plot the value of r in the four scenarios listed at the end of Section 2, for representative values 
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Figure 4: Invariant mass spectra under different scenarios at the two mass values m^: 
(a)m ( - = 170 GeV, (b) = 150 GeV. For each case, the parameters are: solid line: 
scenario (1), sin9 t = 0.0, tan/5 = 1.0, m bL = 300 GeV; dotted line: scenario (2), sin9 t = 
—0.8, tan (3 = 1.0, m bL = 300 GeV; dashed line: scenario (3), sm.6 t = 0.9, tan/5 = 50.0, 
m~ bL = 200 GeV; dot-dashed line: scenario (3), sin# t = 0.0, tan/5 = 1.0, m~ bL = 300 GeV; 
long-dashed line: scenario (4), sin# t = 0.4, tan/5 = 8.0, = 200 GeV. 
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of the parameters, as a function of the stop mass. We see that the value of r is typically 
lower than the top quark value ([13]), that it has a slow dependence on the value of the stop 
mass m, and that it can depend significantly on the stop mixing angle Qf 

The variation of r arises from the competition between the diagrams in Figure [l] in which 
the Goldstino is radiated from the t and b legs and those in which the Goldstino is radiated 
from the W. To understand this, it is useful to think about the limiting cases in which 
each intermediate propagator goes on shell. In the case in which the top quark goes on shell 
in diagrams a,b of Figure [I], the W polarization has the same vlaue ( |13"D as that for top 
decay. In the case in which the b goes on shell, we have the process t — > bW + , for which 
also r = 1/(1 + 2m^/m(). However, the third case in which the x + goes on shell can give 
a very different result. In the limit in which the \ + is pure gaugino, we have the subprocess 
w + —>■ GW + , which leads to purely transversely polarized W bosons. More generally, for 
the process xf ~^ GW + on shell, we have 





V + u\ 


2 + 


V..V2 


2 


2(|V +1 i| 


2 + 


\V- n 


l 2 ) + 


v +12 \ 


2 + \ 




2 



(14) 



where V+, V- are the matrices defined in (0). These individual components vary in impor- 
tance as the masses on the intermediate lines are varied. The role of the chargino diagrams 
in producing a low value of r is shown clearly in Figure [7] Here we plot the value of r as a 
function of the supersymmetry-breaking SU(2) gaugino mass m,2 and observe that r moves 
to a higher asymptotic value as the gaugino is decoupled. 

Beyond this observation, though, the dependence of r on the underlying parameters is 
not simple. As we have seen in the previous section, it is never true that one particular 
subprocess comes almost onto mass shell and dominates the stop decay. This feature of the 
stop decay, which was an advantage in the previous section, here provides a barrier to finding 
quantitative relation between a measured value of r and the underlying parameter set. On 
the other hand, it is interesting that almost every scenario predicts a value of r substantially 
different from the Standard Model value for top decay. 

5 Conclusions 

In this paper, we have discussed the phenomenology of light stop decay through the process 
t — > W + bG. We have shown that this process can be distinguished from t decay through 
the characteristic shape of the bt mass distribution. We have shown also that the fraction of 
longitudinal polarization of the W + in t decays can vary significantly from the prediction ( |i~3"D 
for t. Since these two observables are available at the Tevatron collider, it should be possible 
there to exclude or confirm this unusual scenario for the realization of supersymmetry. 
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Figure 5: Longitudinal W production ratio for the stop t\ decay as a function of mj under 
different scenarios. The dot-dash, dotted, dashed, and dot-dot-dash lines refer, respectively, 
to the chargino scenarios (1), (2), (3), (4) given at the end of Section 2. The two figures 
show (a) tfj-like cases with sin9 t = —0.8, tanf3 = 1.0, and = 300 GeV. (b) pure t^-like 
cases with sin8 t = 0.0, tanf3 = 1.0, and m~ b = 300 GeV. 
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Figure 6: Longitudinal W production ratio for the stop t\ decay as a function of sin9 t under 
different scenarios. The dot-dash, dotted, dashed, and dot-dot-dash lines refer, respectively, 
to the chargino scenarios (1), (2), (3), (4) given at the end of Section 2. The two figures show 
the dependence for (a) m~ tx = 170 GeV, tan(3 = 1.0, and m~ b = 300 GeV. (b) m~ tx = 170 GeV, 
tan(3 = 50.0, and m bL = 300 GeV. 
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Figure 7: Longitudinal W production ratio for the stop t\ decay as a function of the soft 
breaking mass of SU(2) gaugino. The four curves correspond to sin# t = -0.98, -0.8, -0.6, 0.0. 
The other parameters are chosen to be m~ ti = 170 GeV, tan(3 = 1.0, \x = 1000 GeV, and 



m- bL = 300 GeV. 
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